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Abstract 

We consider a mixture of a single-component Bose gas and a two-component 
Fermi gas at temperatures where the Bose gas is almost fully condensed. In 
such a mixture, two fermionic atoms can interact with each other by exchang- 
ing a phonon that propagates through the Bose condensate. We calculate the 
interaction potential due to this mechanism and determine the effective s- 
wave scattering length for two fermions that interact, both directly by the 
interatomic potentials, as well as by the above mentioned exchange mecha- 
nism. We find that the eff'ective scattering length is quite sensitive to changes 
in the condensate density and becomes strongly energy dependent. In addi- 
tion, we consider the mechanical stability of these mixtures and also calculate 
the dispersion and the damping of the various collisionless collective modes 
of the gas. 

I. INTRODUCTION 

The experimental realization of Bose-Einstein condensation in trapped atomic Bose gases 
has stimulated both theoretical as well as experimental efforts to investigate new and 
interesting physics in ultracold atomic gases. In particular, it has spurred a lot of interest 
in achieving a BCS transition in trapped fermionic gases, also because a spin-polarized gas 
of atomic ^Li has theoretically been shown to undergo a transition into a superfluid state at 
a critical temperature that is comparable to those realized in the experiments with trapped 
Bose gases This relatively high critical temperature is due to it's anomalously large 

and negative scattering length of —2160 ao [^]. To reach the temperature regime where 
the quantum degeneracy of the atomic gas becomes important, evaporative cooling is used. 
This technique, however, requires a fast thermalization rate and therefore a high collision 
frequency. For doubly spin-polarized Fermi gases this is not possible due to the exclusion 
principle, and such a high collision frequency must be obtained by using a mixture of either 
Fermi-Bose gases or Fermi- Fermi gases 0. Experimentally, this method of sympathetic 
cooling has already been used to produce the first two-component Bose condensate [Jl0[ and 
has recently also been reported to be succesful for a fermionic two-component mixture 

Mixtures of dilute atomic gases are, however, interesting in their own right, both from an 
experimental as well as from a theoretical viewpoint. Indeed, experimental work on the static 
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and dynamic properties of a binary mixture of Bose gases has included the study of mean- 



field effects [0, relative-phase coherence |]13l, the dynamics of component separation [|14 



and Rabi oscillations . Furthermore, metastable states [^] and quantum tunneling effects 



17| have been observed in spinor Bose-Einstein condensates, and most recently a convenient 



method for the creation of topological excitations in two-component Bose-condensed gases 
has been suggested [|l8l and succesfully carried out experimentally fl^. Theoretical work 
on trapped binary Bose condensed gases, has for example included the stability and static 
properties p0|^2|, the dynamics of the relative phase and the collective modes [Q. Also 
for spinor Bose condensates in optical traps, the phase diagram and collective modes have 



been considered p5|j26[| . In the case of Fermi-Bose mixtures, the density profiles for gases 
confined in a harmonic trap have been studied at nonzero temperatures in a Thomas-Fermi 
approximation [p^ -p9|. In all these cases the theory is based on a mean- field treatment of 
the interactions and neglects the effect of fluctuations. 

Inspired by the well-known physics of ^He-'^He mixtures [^], we here go beyond mean- 
fleld theory and study the effective interaction between two fermions due to density fluctu- 
ations in a Bose condensate. In particular, we calculate the resulting effective interatomic 
s-wave scattering length for two fermions in different hyperflne states with the aim of ma- 
nipulating it in such a way that a BCS transition becomes feasible. The two hyperflne states 
form an effective spin 1/2 system, and therefore the mixture of fermions in two different hy- 
perflne states can be treated as a Fermi gas with spin 1/2. We mostly take the populations 
of the two spin levels to be equal because, if a BCS transition can be achieved at all, this will 
be the optimal situation. Moreover, we also consider the system to be homogeneous, for the 
following reason: If the mixture would be trapped in an isotropic harmonic potential with 
trapping frequency tu, a measure for the overlap of the Bose condensate with the fermionic 
cloud is given by the ratio of the zero-temperature Thomas-Fermi radii for an interacting 
Bose condensed gas and that of an ideal Fermi gas. The former is equal to l{lbNBCiB/^y^^ 



31| , where / = \J?i/muj is the harmonic oscillator length, is the bosonic scattering length, 

and Nb is the number of bosons. The latter equals l{AmFf'^ where Np is the num- 
ber of fermions. For typical experimental parameters and the desirable large numbers of 
bosonic and fermionic atoms, this overlap of the two clouds is rather small. Therefore, in 
order to maximize the effect of the Bose condensate on the interaction between the fermions 
in the mixture, we consider only a spatially homogenous system. Although all the experi- 
ments with Bose condensed gases have been performed with harmonic oscillator traps up to 
now, this is not an unrealistic suggestion since it is certainly possible to create an external 
trapping potential that is more or less a rectangular box . 

Physically, the effect of the Bose condensate on the fermion-fermion interaction is due 
to the exchange of phonons that propagate in the Bose-condensed gas. In order to calculate 
the effect on the interatomic interaction quantitatively, we thus, have to accurately know 
the density-density correlation function in the Bose gas. This is straightforward in the Bo- 
goliubov approximation, whose validity is well-established at such low temperatures that 
the Bose gas is essentially fully condensed. In this manner we obtain the actual interaction 
potential from which we then extract the interatomic scattering length. In addition to the 
effective fermion-fermion interaction, we consider the important question of the stability of 
the three-component system against demixing of the various Bose and Fermi components of 
the gas. We also calculate the excitation spectrum for the Fermi-Bose mixture far below the 
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Bose-Einstein transition temperature, where the dynamics of the gas is in the colhsionless 
regime. Theoretically, this amounts to doing a so-called RPA (Random Phase Approxima- 
tion) calculation for this mixture. As a result we find not only the eigenmodes of the gas, 
but also the Landau damping of these modes due to the imaginary part of the 'RPA-bubble' 
diagram. Note that from a fundamental point of view the stability and the excitation spec- 
trum are strongly related, because a signature of the demixing-inst ability is the occurrence 
of a mode with a purely imaginary frequency. 

The paper is organized as follows. In Sec. |I| we first calculate the effective interaction 
and scattering length of two fermions in the presence of a Bose condensate. In particular, 
we show that for realistic conditions, the scattering length strongly depends on the collision 
energy of the atoms, which is important when considering the prospects of a BCS transition 
in a spin-polarized potassium gas. In Sec. |T| we then consider the stability of the three- 
component system and show that in general this does not lead to very stringent constraints 
on the densities. Finally, in Sec. |V], we calculate the long-wavelength collective mode 
spectrum and damping for the Fermi-Bose mixture. We end with a summary in Sec. 0. 



II. EFFECTIVE INTERACTION 

In this section, we calculate the effective interaction and scattering length of two fermions. 



The calculation is performed by means of functional methods [Q, because even in the 



presence of a Bose condensate, the RPA calculation for the collective modes can then be 



performed in a purely algebraic manner in Sec. IV and avoids the complications of an explicit 



evaluation of Feynman diagrams. Moreover, it gives some more insight in the physical nature 
of the colhsionless collective modes. Of course, it is also possible to do the same calculation 
in the operator formalism. If performed correctly, it gives identical results. 



A. Theory 

In accordance with the previous remarks, we thus start from the functional-integral 
expression for the grand-canonical partition function of the mixture. It reads 

Zg, = I f/[0*]rf[0]f/[v^*]c?[v^]exp|-i(5B[0*,0] + ^^[V'*,V'] + 5/[r,0,v^*,v^])} , (1) 

and consists of an integration over a complex field 0(x, r), which is periodic on the imaginary- 
time interval [0, hjS], and over the Grassmann fields ipai'^, t), which are antiperiodic on this 
interval. Therefore, 0(x, r) describes the Bose component of the mixture, whereas ^/'q,(x, r) 
is associated with the Fermi components. For the latter we actually need two fields, because 
the fermionic atoms can be in either one of the two hyperfine states Furthermore, the 
total action of the mixture consists of a term for the Bose-gas 




a term for the Fermi-gas that accounts for the fact that the Pauli principle forbids s-wave 
scattering between fermionic atoms in the same hyperfine state 
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+ ^|^„(x,r)n^_„(x,r)r| , (3) 

and a term describing the interaction between the three components of the Fermi-Bose 
mixture 

Si[<P*,<P,r,i^]=T.Ta / dr / (ix|^,(x,T)n0(x,T)|2 . (4) 

In these expressions we have introduced the two-body boson-boson T(ransition) matrix ele- 
ment Tb — 4:7rh'^aB/mB, the two-body fermion-fermion T-matrix element Tp — Ai^Ti^aF /mp-, 
and the two two-body T-matrix elements = 21^11? aa/ruR that describe the interactions 
between a boson and a fermion in the spinstatc \a). Here a = {], 1} denotes the hypcrfine 
components of the Fermi-gas in the effective spin language. In addition, fi denotes the chem- 
ical potential of the Bose-gas and /la denotes the chemical potentials of the two components 
of the Fermi gas. Note that the latter do not need to be identical, because the Fermi gas 
is in general not in equilibrium in spin space, due to the generally slow relaxation rates 
between the hyperfine degrees of freedom. Finally, the masses of the bosonic atoms and the 
fermionic atoms arc denoted by ms and mp, respectively. The reduced mass is denoted by 
f^R = mpmB/ {mp + m^)- 

Since we only consider the gas at such low temperatures that the Bose gas is essentially 
fully Bose condensed, we proceed by performing the usual Bogoliubov substitution for the 
Bose-fields, i.e., (j) = ^/nB + (f)', and neglect all terms of higher than second order in 0' or 
0'*. Note that we have implicitly also neglected the depletion of the condensate and put the 
condensate density equal to the total density of the Bose gas. For a weakly- interacting 



gas with ynBa% <^ 1 and at low temperatures, this is clearly justified. In this manner we 
end up with the Bogoliubov approximation Sb[4>*:4'] — hpVTBnB/'i + Sb[4>'*: 4>']: where V 
is the total volume of the gas and the action for the fluctuations can be written as 



SbW*,(P'] = --J^ dTdr'J d^d^ct>'\-^,T) ■ G-^(x,T;x',r') • </>'(x',t') , (5) 



2 Jo 

if we introduce the vector field 



"(x.t) 



and the corresponding Green's function G(x, t;x',t') that obeys 
— ^G~^(x, r; x', r') 

( hdr - h'^V^/2m + TBnB Tub \ a .-x 

Tub -^a.-^W2m + T^nJ'^(^-^)'^(^-^)- 

Note that the linear terms in (f)' and (j)'* have dropped out of the action because of the 
Hugenholtz-Pines relation 
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/i = TbUb + TaUa , (8) 
a 

which also incorporates the mean-field effects due to the nonzero spin densities in the 
Fermi gas. The same relation is also used to eliminate the chemical potential from Eq. (|^. In 
addition, the interaction term in the total action becomes in the Bogoliubov approximation 

with a 'current source' defined by 

J{^,T) = ^YTa\M^,r)\'(\) ■ (10) 



To include the effect of phonon exchange on the interaction between two fermions, we 
now need to integrate over the Bose fields. This is straightforward, since it only involves the 
evaluation of a gaussian integral. The result for the total effective fermion action is thus 

S'^^hp*, i;] = Splip*, ^] + — drdr' / (ixrfx'J^(x, r) ■ G(x, r; x', r') • J(x', t') 

2h Jo J 

= SF[r,^P] + 2^7 drdr' J dxrfx'|^„(x,r)|V,,„,(x,r;x',r')|^„Kx',r')P , (11) 

where VQ,a'(x, t; x', r') is the effective interatomic potential due to phonon exchange. At 
this point, it is important to note that to be able to use in Eqs. (0), and (|) the two- 
body T- matrix elements, instead of the real interatomic (singlet/triplet) potentials, we have 
already integrated out fluctuations with momenta higher than a certain cutoff hA. Thus, 
in principle there is a cutoff on all the momentum integration in the rest of this paper. We 
will comment on the effects of this shortly. Hence, inverting Eq. (|^) by means of a Fourier 
transformation, we find that ^^^^/(x, r; x', r') is given by 



K,a'(x,r;x',r') 

6(k) 



-2TaTa'nBYl / 



<A (27r)3 



^|[^gik-(x-x')-i<^„(T-T') 



{huJnY + e{}i){eiy) + 2TBnB) 



'12) 



where e(k) = ^^k^/2m and uon = 2TTn/hp are the even Matsubara frequencies that ac- 
count for the periodicity of the Bose field 0(x, r) and therefore of the Green's function 
G(x, r; x', r'). The phonon-exchange mechanism thus also induces an interaction between 
fermions with the same spins. Due to the Pauli-principle, however, this interaction can again 
at best be of a p-wave nature, and is in general negligible. As a result, we from now on 
only consider the contributions to the interaction potential between particles with opposite 
spin. We also neglect the frequency dependence in Eq. (0), and consider only the static 
contribution. This implies that the relevant collision energies of the fermions must be much 
less than the bosonic mean-field interaction. 

The resulting instantaneous potential is given by \^(x, r; x', r') ~ V^^dx — x'|)5(r — r') 
with 
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Here we have defined the coherence length ^ = 1/ ^/XQwubc^b in the Bose-condensed gas. The 
effective potential depends on three parameters, the cutoff A, the coherence length ^ and 
uiBa^aJm^aB- As we will show now, however, there is a separation of length scales, and 
1/A is in general much smaller than the other relevant length scale, the coherence length ^. 
This means that we can safely ignore the cutoff and take the limit A — >■ cxo. To estimate the 
magnitude of the cutoff, we make use of the fact that it is determined by the requirement that 
the bare interatomic interaction has renormalized to the two-body scattering matrix. This 
renormalization of the bare interaction is described by the Lippmann-Schwinger equation 

Ta ~ K) ^ A<fc<Aa (27r)3 2e(k) ' ^ ' 

Here Aa denotes the ultra-violet cutoff provided by the interatomic potential, and A is 
the momentum scale up to which fluctuations in the Bose gas are integrated over. The 
bare interaction Vq is chosen such, that in the limit A = the result does not depend on 
the high-momentum cutoff A^ [^. Requiring the renormalized interaction to be within 
approximately a fraction x from the two-body transition matrix leads to 

A ~ ^ . (15) 

If we compare this value for the cutoff with the coherence length ^, the product of the two 
is given by 

eA . (16) 

SJiinBa% 



Due to the presence of the factor \JnBa% in the denominator, this quantity is in general 
much larger than 1, and we can safely take the limit A ^ oo. The effective interaction in 
this case becomes simply 



fi^mBCi^a 
4m%aB^^ 

and has the form of a purely attractive Yukawa potential 



V^^ir) = _l_!!lE^,-r/, ^ (17) 



B. Results 

With this interaction potential V^^(r), we want to determine the effective s-wave scatter- 
ing length a^p for two fermions interacting also through the mechanism of phonon exchange. 
We can associate an s-wave scattering length a'^p with the effective interaction V'^^{r) as 
follows 

af = -\iJ-^. (18) 

k[0 k 
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Here So{k) denotes the phase shift of the partial wave with angular momentum / = 0. The 
phase shift 5o(^) is defined in terms of the asymptotic form 



lim Uo{r; k) sin(A;r + So{k)) (19) 

for the / = partial wave uo{r;k), that can be calculated from the radial Schrodinger 
equation 

^ +^V^''{r)-kAuo{r;k) = 0. (20) 



At this point, we have to keep in mind that there is already a scattering length ap for the 
fermions due to the interatomic potential. We can take this interatomic scattering length 
into account by imposing suitable boundary conditions on the partial wave uo{r;k), such 
that we recover ap H the interaction due to phonon exchange vanishes, i.e., if V^^{r) = 0. 
The boundary conditions imposed are such that the derivative and the magnitude of the 
wave function at r = are equal to those of the function sin[A;(r — ap)]. In the case of a 
positive scattering length, an alternative procedure would be to add a hard-core potential 
with a range to V^^i^r). Both these procedures are justified because the range of the 
interatomic potential is much smaller than the range of the effective potential, and we have 
checked that numerically they indeed yield the same results. 

The phase shifts due to the effective potential V^^{r) for mixtures of ^"K and ^Li with 
^''Rb, are show in Figs. |I| and H, respectively. The interatomic scattering lengths for ^°K-^'^K, 
^Li-^Li and ^''Rb-^^Rb collisions are taken to be 160 ao, —2160 ao and 109 ao, respectively 
@;36]. The scattering lengths for and ^Li with *^Rb are, as far as we know, not known 



and we have taken them to be equal to 100 ao, which is presumably a typical value. In Fig. |^ 
also the phase shifts for a mixture of "^^K and ^^K are shown, where the scattering length for 
and ^"^K-^^K collisions has been taken to be 5 ao and 1000 ao, respectively The 



various lines corresponds to different values of the condensate density. It is clear from these 
figures that for a given mixture, the phase shifts as a function of momentum asymptotically 
all have the same slope independent of the value of the condensate density. This slope 
corresponds to the interatomic scattering length ap. However, at long wavelengths the phase 
shifts can be significantly different from the asymptotic limit, and can depend strongly on 
the coUisional energy. For a mixture of ^'^K-^''Rb this is certainly the case. The effective 
interaction can even become attractive, instead of repulsive, around the Fermi momentum, 
which opens the possibility for a BCS transition to a superfluid phase just like the exchange 
of phonons in the lattice leads to superconductivity in metals For a mixture of ^Li- 
^''Rb the effects are less pronounced, and amount to a slight enhancement of the already 
very large and negative scattering length. The case of "^^K-^^K may be exciting, because of 
the possibility of a rather large scattering length for ^"K-^^K collisions, which can be of the 
order of 1000 ao If this value is correct, there are resonances in the effective scattering 



length as shown in Fig. ^. This strongly resembles the existence of Feshbach resonances in 
the scattering length as a function of the applied bias magnetic field |^ . 
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III. STABILITY 



To be able to interact, the Bose and the two-component Fermi gases have to be overlap- 
ping. Therefore, we want to consider the stability of the gas against demixing of the Fermi 
and Bose components. To determine the stability of the mixture, we consider the matrix of 
second-order derivatives with respect to the densities of the free energy F = Fp + Fb + Fj, 
which consists of a fcrmion term Fp. a boson term Fb and an interaction term Fj. For the 
mixture to be stable we have to require that all eigenvalues of this matrix are larger than 
zero, i.e., that the free-energy surface is convex. At the low temperatures of interest to us, 
the free-energy density is equal to the average energy density (E) /V and we get 



{E)/V = {Ef)/V + {Eb)/V+{E, 



)/V 



^(Gn'mnf + nf)-^ + Tfn.n, + Tb^ + {T,n, + T,n,)nB 



(21) 



Here, the first term is the kinetic energy of the two-component Fermi gas, the second and 
third terms are the interaction energies of the individual Fermi and the Bose gases, respec- 
tively, and the fourth term is due to the interaction between these gases. We have neglected 
the kinetic energy of the Bose gas because of the low temperatures and densities. Note that 
the effective fermion-fermion T- matrix element T|? = Amh^a^p jmp has to be used, that in- 
cludes the effect of phonon exchange. This is the case because, as long as the Fermi energy 
is much less than the average mean-field energy of the Bose condensate, we have to include 
the effect of fluctuations, which change the static properties of the mixture by renormalizing 
the fermion-fermion interaction as we have just discussed. 
The matrix of second-order derivatives is given by 



d^F 

d^F 



d'^F 



dn-rdrtB 
d^F 

dn^drtB 
92 F 



duBdny dnsdnB / 



Assuming that we start in the stable part of the phase diagram, the onset of an instability 
is signaled by the point where its determinant becomes equal to zero. Therefore, requiring 
the determinant of this matrix to be larger than zero is a sufficient condition for stability. 
It reads 



TT 



TT 



rrLBmp 



l^F , 



m 



R 



TT 



/ci a' 



at + 



TT 



„eff 



+ 8 



a|a| rriBmp 



- 16 > 



(22) 



R 



where = (GTr^nQ)-*^/^ denotes the Fermi momentum associated with the spin state \a). 
Note that a^p is a function of the condensate density hb- In evaluating the derivatives of 
the free energy with respect to the density of the Bose gas, we have not taken this implicite 
dependence into account, which is sufficiently accurate for our purposes as long as we are 
not too close to a resonance. Near a resonance a more involved treatment is necessary, also 
because the phaseshift is then very strongly momentum dependent. 
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The surface in the ria — ub volume where the equahty sign holds is called the spinodal 
surface. The spinodal surface divides the phase space into a region where the mixture is 
(meta)stable, and one where it is unstable and separates into distinct phases in the stable 
part of the phase space. If the effect of phonon exchange is optimized by putting = 
we have ki = = kp, and Eq. ([2^ ) becomes 




(a? + (^i. 


^ msmp 


( 7T 






{kpaf 


(a| - ai 




( ' 









+ 2 > . (23) 



If we also put the condition further simplyfies to 



TT 




kpa^p J \kpa'^p aBd^p / 



(24) 



The first factor in the left-hand side of Eq. (^) corresponds to the demixing of the two 
fermion components of the gas 0] , whereas the second factor represents the demixing of the 



Bose- Einstein condensate and the fermion components |]39|. Note that it depends on the 
various masses and scattering lengths involved, which of the two instabilities occurs first. 
For equal scattering lengths, i.e., a'^p = as = = ai = a, the demixing of the Bose and 



Fermi gases always occurs first, and we reproduce the result of van Leeuwen and Cohen 



vr niB mp 
> — + — . 25 

kpa nip TUB 

In our numerical calculations, the scattering lengths and a| have always been taken 
equal to each other and Eq. ( plf ) applies. In the case of a mixture of and ®^Rb this 
condition roughly leads for the total fermion density np only to the restriction np < 10^^ 
cm~^. For a mixture of ^Li and ^''Rb the condition on the total density becomes np < 10^^ 
cm~^, and for a mixture of and ^^K it reads np < 10^° cm~'^. The latter condition seems 
to be quite restrictive. However, if we take instead of ub — 5 ag a different value that is 
within the present uncertainty for this scattering length, i.e., as = 25 ao, the condition for a 
mixture of and '^^K becomes only np < 10^^ cm~^, which is much more favorable when 
one considers the prospects of achieving a BCS transition in this case. The reason is that in 
BCS theory, the critical temperature is given by Tc = {Sepe'^'^ /kBn) exp {—TTcot{6o{kp)) / 2} , 
with 7 Euler's constant. From this expression it follows that if the densities are low, the 
critical temperature for the BCS transition is also very low. Therefore, if the results of 
Sec. H for a and ^^K mixture are to be of use, as far as achieving a BCS transition is 
concerned, it is crucial that relatively high densities are realizable. 



IV. DISPERSION OF COLLECTIVE EXCITATIONS 

We next want to consider the excitation spectrum of the gas. The collective excitations of 
the mixed gas are coupled modes of the fermionic spin densities, and the Bose condensate. It 
is therefore convenient to perform a Hubbard-Stratonovich transformation to these fermionic 
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densities This amounts to introducing two real auxiliary fields P|(x, r) and pj,(x, r), by 
rewriting in the integrant of the functional integral for the partition function the factor 
due to the fermion-fermion interaction as, 

expj-^^ dT J rfx|^/'|(x,r)|Vi(x,r)p| = J d[p^]d[pi\ 

For reasons that become clear shortly, we denote the resulting action for the fermions by 
Sp[%lj*,ip\. It is quadratic in the fermion fields, and reads 



S^[rM = llfJdrjdJra{^,T 



^T. 7^ -pc. + Ti.p_«(x,r) 

OT imp 



+ Tc,|0(x,r) 



^a(x,r) 



(27) 



Upon integrating over the fermionic fields, the grand-canonical partition function describing 
the gas is equal to a functional integral over the Bose field and the density fields only, with 
an effective action that reads 



a 

Jo J ^^^^^T(x.^)Pi(X'^) ' 



(28) 



5^«h,Pi,r,0] = 

in terms of the Green's functions 
G;i(x,r;x',r') 

= - ^ - P. + T^P-.(x, r) + Tj0(x, r)p| <5(x - x')5(r - r') . (29) 

Expanding this action around its minimum, by requiring the linear terms in the fluctua- 
tions to be zero, will result in the Hartree approximation for the equilibrium densities. This 
is sufficiently accurate because we are dealing with a fermionic gas in a nonmagnetic phase, 
where the Fock term in the selfenergy is zero, due to the spin-symmetry of the action. It 
also explains the use of the symbol H(artree) in Eq. (|27|). Inserting thus into the right-hand 
side of Eq. (53) 



p„(x,r) = + p^(x,r) 



and 



x,r 



(30) 



(31) 



we obtain from the requirement that the linear terms in the fluctuations vanish again the 
Hugenholtz-Pines relation 
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- /i + TbUb + TaG^ (x, r; X, r) = , (32) 

a 

and in addition the expected equation for the average spin densities 

= (x, r; X, r) , (33) 

where 

^^"(x, r; x', r') = |/^|: - ^ - + T^P-a + T^ub^ 5(x - x')5(r - r') (34) 

is the usual fermionic one-particle propagator, or two-point correlation function, in the 
Hartree approximation. We, therefore, recognise in Eqs. (|32|) , (|33|) , and (Q) the selfconsis- 



tent Hartree equations for the boson and fermion densities at given chemical potentials /x, 
and 

To find the theory describing the fluctuations around this equilibrium, we have to perform 
the expansion around this minimum up to second order in the fluctuations p'^ and 0'. The 
poles in the Green's function of the resulting theory give us the desired dispersion of the 
collective modes, for they are also the poles in the linear response of the densities to an 
external perturbation. If we introduce again the vector notation 



A,. ^ I ,f I , (35) 



and also 



P'^n = ^J^*^'" , (36) 

the quadratic part of the effective action for the fluctuations and can be conveniently 
written in momentum space as 

c^(2)r.A' n'l = i { '^M V i M<^<^(k,«^n) M<^p(k,2u;„) \ / „ \ , . 

where we have defined the matrices 

, , _ / -ihun + e(k) + E(k, zw„) E(k, iw„) \ , . 

lUn) - I ^^^^ ^^^^ ^^^^ ^ ^^^^ ^ ^^^^ ^^^^ I , {6^} 



and 



and 



M„(M.„)^T. (39) 



M^^(k, lu;^) - M^,(k, lu^r.) - TbVt^ ^.^^^^^^ ^^^^ ^^^^^^^ ^^^^ j • (40) 
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We have also introduced the selfenergy of the Bose gas 



E(k, iun) = TbUb + T?ni-(k, iun) + Tpn|(k, iujn 



riB , 



(41) 



and the RPA-bubble, or equivalently, the density-density response function of the ideal Fermi 

gas 



n„(k,iu;„) 



dp n(ea(k + p)) - n(eQ(p)) 



(42) 



where the dispersions obeyeQ,(k) = h^k'^ /2m+Tpp_a+TanB—fJ,a, and ?T,(e) = l/(exp(/5e) + l) 
is the usual Fermi-Dirac distribution function. Note that S(k, ia;^) and Ila(k., iuJn) are 
invariant under the substitution (k, icj^) (— k, —iun), due to the time-reversal symmetry 
of the problem. 

In the long- wavelength limit, i.e., for small k, Eq. (|4^) can be rewritten as 



n„(k,iu;„) 

where 



mka 



de 



1 — arctan 



^l3{tc,e-fJ,ci+Tpp-a+TanB) 



+ 1 



(43) 



(44) 



and we have defined the Fermi energies and wavevectors by ea = Pa — Tpp-a — T^ns = 
h'^k'^/2m. The analytic continuation of Eq. (|43| ) to physical energies lu = iuJn reads 



2h^Ti'^ Jo dt 



~ , dur 
ae- 



1 log 

2 ^ 



vr 



(45) 



with Xa = muj/hkak. We want to find the zero-temperature result and the lowest-order 
corrections in the temperature. This can be done by means of a Sommerfeld expansion, 
which amounts to expanding the real part of the expression between square brackets in 
Eq. ( ^51) around e = 1. Doing so, we find in first instance 



n„(k,u;) 



de- 



2h\^ Jo de 

^ai^a ~ 3) 



log 



1 



:i-x2)2 



2 



X ~\~ 1 



Xn 



Integrating then over e, we obtain for the real part in lowest order 



nl°)fk. 



UJ 



mkana{0) 
2h\^ 



Xr 



1-ylog 



Xa + l 



Xqi 1 



and for the imaginary part exactly 

Im[n.(k,^)] = "2Sf^2M0) '^"'""('^-'') • 



(46) 



(47) 



(48) 
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At zero temperature Eqs. 
temperature RPA bubble [42 



and (^81) reduce to the well-known result for the zero- 
Note that the imaginairy part in Eq. (|48|) is valid for all 
temperatures and cannot simply be expanded as a power series in ksT /ep. This is not true 
for the real part of n„(k, cj) and the lowest-order temperature correction is given by 



3) 



(49) 



To find the long-wavelength dispersion, i.e., u = ck or equivalently Xa = mc/hka, of the 
collective modes of the mixture, we require that the determinant of the fluctuation matrix in 
Eq. (|37D equals zero. Calculating the determinant, the terms of 0{k^) drop out, indicating 
that the collective excitations are indeed gapless as assumed by our ansatz u = ck. Leaving 
out terms that are of O^k"^), because they do not affect the linear part of the dispersion 
relation, and dividing out an overall factor of k"^, we ultimately find the result 



T|n|(c)n|(c) - l] [-mc^ + S(c 
-TpriB [TfU^{c)T^U^{c) - T^n^(c)] T^U^{c) 



. 



(50) 



Here, we have for convenience introduced the shorthand notation liaic) = na(k, c/c). If we 
put T| = Tj = 0, this simply gives 



-mc 



S(c) 







(51) 



The first factor describes the collective modes of the Fermi gas and the second factor de- 
scribes the Bogoliubov modes of the condesate, which are of course decoupled in this case. 
For small values of T| and T| these modes also exist, but the dispersion is changed. The 
lowest-order correction in the boson-fermion scattering lengths to the Bogoliubov speed 
of sound Co = jTsnB/m is determined by 



mci 



S(co) 



TfHb [Tj.nT(co)T;n^(co)-TTnt(co)]T;n^(co) 

inMconiico) - 1] 
TpUB [TFn|(co)T|n|(co) - T|n^(co)] T|n|(co) 



(52) 



To find the full solution we can just iterate Eq. ( p2D and the result converges rapidly to a 
solution of Eq. (|50|) . Note that Eq. (^) has an imaginary part and therefore also describes 
the damping of the Bogoliubov mode. Under the circumstances that we have studied in 
Sees. H and |T| the correction on the speed of sound due to the presence of the Fermi 
gas is small and at most about 5% of the uncoupled Bogoliubov result. This is important 
for our purposes, because it shows that for the calculation of the effective fermion-fermion 
interaction, we do not need to selfconsistently include the effect of the Fermi gas on the 
density fluctuations of the condensate. 

To find the other propagating solutions of Eq. (^) we need to be more careful. For 



clarity, we first again treat the case where T| 
solutions of 



T| = 0. The zero-sound modes are now the 
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[cmic) - 1 







(53) 



The product n|(c)n|(c) diverges logarithmically for two value of c, i.e., for the Fermi veloc- 
ities Cq = hka/m. In principle there are therefore four solutions to this equation, which can 
be found by expanding around either or C|. Without loss of generality, we assume that 
C| > Cj. To find then the solution that is not overdamped, we have to expand around C| and 
find the solution which is slightly bigger than c-f, i.e., c = + 6c with 6c > 0. In this way 
we make sure that the imaginary parts of both Il^{c) and nj(c) are equal to zero at zero 
temperature. Experimentally we are always in the weak-coupling limit, which implies that 
Tp <^ h'^ /mk^. In this case the zero-temperature expressions for n|(c) and n|(c) become 
explicitly 



n°(ct + 6 c) 



mkf 



1 ^ 



+ 0{6c) 



(54) 



and 



n;(ct + 6 c) 



mki 



^) + l 



+ 0{6c) 



(55) 



respectively. In the limit that (c| — c^) ^ 5c we thus get 



6c = 2c| exp 



mkfTp 



(56) 



where we have made use of the fact that in this limit kj^ — k-^ = 0{6c). Note that our result 



differs from that of Fetter and Walecka ||4^. They use an approximation that does not obey 
the Pauli exclusion principle, because there is s-wave scattering between particles that are 
in the same spin-state. As a result, their mean-field energy for an atom in the spin state |a) 
is TpJ2a^a instead of Tpn_a. For = this effectively means a factor of two reduction 
of the interaction strength. In the limit that (c| — Cj) ^ 5c we find at zero temperature 



6c = 2ct exp 



47r2^2 



mk^T^U^ic^) 



-2 



(57) 



Let us now consider the effect of the presence of the Bose condensate and do the same 
calculation for the mixture of the Bose condensed gas with the two-component fermion gas. 
In the limit that (c| — Cj) ^ 5c we again find 



6c = 2c| exp 



2t2 



mkfTp 



(58) 



This is the same result as in Eq. (|56| ) where the fermionic and the Bose condensed gas are 
decoupled, and is due to the fact that in this limit the zero-sound mode is a pure spin wave. 
The density profile is therefore constant and, since the fermionic and the Bose condensed 
gas couple only through the density fluctuations, there is no effect of the presence of the 
Bose gas. On the other hand, in the limit that (c| — C|) ^ 5c we now obtain 
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2*2 



5c = 2c| exp 



where we have defined the quantity 

Tf 



mktTpX 



(59) 



X 



+ TBnB + TfUl{c^)nB] 



r2n0(ct)(-mc2 + TbUb) - T^riB + 2TFT^T^Yl\{c^)nB 



(60) 



Hence, Eq. (|59|) reduces to Eq. (|57D if we take = T| = 0, as it should. The same is of 
course true if = 0. Note that Eq. (|^) offers the opportunity to obtain also the small 
nonzero-temperature corrections to the above results, if this turns out to be necessary for a 
particular application. 



V. CONCLUSIONS 



We have calculated the effect of phonon exchange on the scattering length for two 
fermions with different spins. Under the right circumstances this effect can be quite large 
and is then a possible way to experimentally tune the fermion-fermion scattering length. 
This may in particular be useful for the achievement of a BCS transition in a mixture of 
and ^^K, although it appears that a more precise determination of the various scattering 
lengths involved in this case is necessary to make sure of this. In addition, we have analyzed 
the stability and the mode structure of the Bose-Fermi mixture. Our results in the latter 
case, which are valid for an arbitrary ratio of the densities of the three components, agree 
in limiting cases with expressions obtained by other authors E3 . 
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FIGURES 



FIG. 1. The phase shift as a function of the momentum k, for a mixture of K and Rb. 
The condensate density is (1) = 1 x 10^° cm~^, (2) = 1 x 10^^ cm~^, and (3) = 1 x 10^^ 
cm~^. 

FIG. 2. The phase shift Sq as a function of the momentum k, for a mixture of ^Li and ^^Rb. 
The condensate density is (1) = 1 x 10^° cm"^, (2) = 1 x 10^^ cm-^, and (3) ns = 1 x 10^^ 
cm~^. 

FIG. 3. The phase shift Sq as a function of the momentum k, for a mixture of and ^^K. 
The condensate density is (1) ns = 4 x 10^^ cm-^, (2) = 5 x 10^^ cm'^, (3) „^ = g x 10^^ 
cm-^, (4) ns = 7 X 10^^ cm-^, (5) ub = 8 x 10^^ cm-^. 



FIG. 4. The scattering length ofjp clS 8; function of the density of condensed atoms ub, for a 
mixture of and ^^K. 



18 



M.J.BijIsma etal. fig.1 



0.02 - 



0.00 



-0.02 - 



-0.04 - 



-0.06 




0.00000 



0.00010 0.00jD20 

k (in units of ) 



M.J.BijIsma etal. fig. 2 



1.50 




M.J.BijIsma etal. fig.3 



2.0 




M.J.BijIsma et al. fig.4 




